Milstein's type schemes for fractional SDEs 



Mihai GradinarrQ and Ivan NourdirH 

Resume. On ctudie la vitesse exacte de convergence de certains schemas d'approximation associes a 
des equations differentiellcs stochastiques scalaircs dirigccs par le mouvcmcnt brownicn fractionnaire 
B. On utilise le comportement asymptotique des variations a poids de B, et la limite de l'erreur entre 
la solution et son approximation est calculee de facon explicite. 

Abstract: Weighted power variations of fractional Brownian motion B are used to compute the 
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exact rate of convergence of some approximating schemes associated to one-dimensional stochastic 
differential equations (SDEs) driven by B. The limit of the error between the exact solution and the 
considered scheme is computed explicitly. 
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Ph ; 1 Introduction 
-C 

Let B = (Bt)te[o,i] be a fractional Brownian motion with Hurst index H E (0, 1). That is, B 
is a centered Gaussian process with covariance function given by 

Cov(B s ,B t ) = \(s 2H + t 2H -\t-s\ 2H ), a ,t€[0,l]. 

> 

For H = 1/2, B is a standard Brownian motion, while for H 7^ 1/2, it is neither a semimartin- 
', gale, nor a Markov process. Moreover, it holds, for any p > 1: 

O" E\B t - B s \ p = c p \t - s\ pH , s,t€ [0,1], with c p = E(\G\ p ), G ~ ^(0, 1) 

and, consequently, almost all sample paths of B are Holder continuous of any order a E (0, H). 

The study of stochastic differential equations driven by B has been considered by using 
several methods. For instance, in [22] one uses fractional calculus of same type as in [25]; in 
[2] one uses rough paths theory introduced in [TT], and in [19] one uses regularization method 
used firstly in [2"5] . 

In the present paper, we consider the easiest stochastic differential equation involving 
fractional Brownian motion, that is 



X 



dX t = a(X t )dB t , t G [0, 1], X = x£R. (1.1) 

Here and in the rest of the paper, a G C°°(M) stands for a real function which is bounded 
with bounded derivatives. Let us denote by (j> : M 2 — > R the flow associated to a, that is the 
unique solution to 

rv 

4>(x,y)=x+ / a((f>(x, z))dz, x,yEBL (1-2) 
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Assume that the integral with respect to B we consider in (|1.1[) verifies the following Ito- 
Stratonovich type formula: 

f(B t ) = /(0) + f f'{B s )dB s , [0,1], / : R -» R smooth enough. (1.3) 

Jo 

Then, combined with (II. 2p . one easily checks that 

X? = </>(x,B t ), t€[0,l], (1-4) 

is a solution to (jl.ip . 

Approximating schemes for stochastic differential equations of the type 

dX t = a(X t )dB t + b(X t )dt, t E [0,1], X = ieR, (1.5) 

have been considered only in few articles. The first work in that direction is [10]. Precisely, 
whenever H > 1/2, it is shown that the Euler approximation of equation (|1.5p - but in the 
particular case where cr(X t ) is replaced by cr(t), that is the so-called additive case - converges 
uniformly in probability. In [18] one introduces (see also [21] ) some approximating schemes for 
the analogue of (jl.5p where B is replaced by any Holder continuous function. One determines 
upper error bounds and, in particular, these results apply almost surely when the driving 
Holder continuous function is a single path of the fractional Brownian motion B, and this for 
any Hurst index H E (0, 1). In [12], upper error bounds for Euler approximations of solutions 
of (II. 5p are derived whenever H is bigger than 1/2. The convergence of Euler schemes has 
also been studied in [3] in the context of the rough paths theory. 

Results on lower error bounds are available only since very recently: see [13] for the 
additive case, and [15] for equation (jl.5p (see also [H] where approximation methods with 
respect to a mean square error are analysed). More precisely, it is proved in [15] that the 
Euler scheme X = {X^} n ^ associated to (|1 . 5j) verifies, under some classical assumptions 
on a and b and whenever H E ( l /2, 1), that 



n 2H ~ l 



X[ n) - ZiJ - ^ / a'(X s )D s X 1 ds, as n 00. (1.6) 
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Here, D S X\ denotes the Malliavin derivatives of X\ with respect to B. Observe that the upper 

and lower error bounds are obtained from an almost sure convergence, which is somewhat 

surprising when compared with the case H = 1/2, see below. In |15| . it is proved that, for the 

i n \ 

so-called Crank-Nicholson scheme X = {X } nG N associated to (jl.ip and defined by 

[ X< (e+i)/ n = X £/l + \ (°"( X l/n) + cr ( X (£+l)/n)) ( B (i+l)/n ~ B £/n), ^ E {0, . . . , n - 1}, 

(1.7) 

the following convergence holds for a regular enough and whenever H E (Y3, Y2): 



n a 



x^-Xi 



^5o, Va < 3^ - 1/2, asn^oo. (1. 



In the particular case where the diffusion coefficient a verifies cr(x) 2 = ax 2 + /3x + 7, for some 
a,P, 7 E M, one can derive the exact rate of convergence and one proves that, as n — > 00: 



Xi n) - X\ 



l ^^a{X 1 )G. (1.9) 
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Here, G is a centered Gaussian random variable independent of Xi, whose variance depends 
uniquely on H. In particular, the upper and lower error bounds are obtained here from a 
convergence in law. 

As we said, the convergence in (|1.6p is somewhat surprising, since there is no analogue for 
the case of the standard Brownian motion. More precisely, when H = l /2, it is proved in [8j 
that the Euler scheme (jl.6p verifies (by denoting X lto the solution of (jl.5p in the ltd sense), 

If 1 

V^x[ n) -X{ td ^-J-Yi a{X l * )a'(Xf°)Y- x dW s , as n -> oo. (1.10) 
L J v2 Jo 

Here, W is a Brownian motion independent of the Brownian motion B and 

Y t = exp (J b'(Xl to )ds - ±J a</(X s Ito ) + a'{Xf°) dB^j , t G [0, 1]. 

On the other hand, it can be proved (see Remark 14.21 2) that, for the Crank-Nicholson scheme 
p.7j) . we have, as n — > oo: 



n 



Law 1 



1 lja 2 )"(X^)[VT5dW s + 3dB s } + ^ 1^ a(a 2 )'"(Xf*)ds, (1.11) 



where X Str denotes the solution of (jl.ip in the Stratonovich sense. 

In the present paper, we are interested in a better understanding of the phenomenoms 
observed in (|1.9|) . (jl. 1(J|) or (|1.11|) . What type of convergence allows to derive the upper 

and lower error bounds for some natural scheme of Milstein's type? More precisely, let us 
define, by induction, the family of differential operators (S^ ).,eNu{o} as 

®°f = f, ® X f = f<J and, fori > 2, & f = ®\&~ l f). (1.12) 
For instance, the first @ J cr's are given by: 

9°<r = a, 9 x a = at/, 9 2 a = aa' 2 + a 2 a", 9 3 a = aa' 3 + 4a 2 a'a" + a 3 a"', etc. 
Now, let us consider the following scheme introduced in [18J: 



x<g = x^ + E7 =0 jjhY. ( A %n?' +1 , l G {0, • • • , n - 1}, 



;i.i3) 



the integer m G N U {0} being called the size of X = {X^} n& ^. Here, for j, n G N and 
£ G {0, . . . , n — 1}, we set ABe/ n instead of B(i+i)i n —Bu n for simplicity. The idea of introducing 
these schemes will be better explained in Section 3 below. For the moment, just observe that 
Euler (resp. Milstein) scheme corresponds to m = (resp. m = 1). 

The aim of the present paper is to answer the following questions. Does the sequence 
{X^} nS N converge? Is the limit Xf given by (|1.4|) . as could be reasonably expected? What 
is the rate of convergence? Are upper and lower error bounds obtained from a convergence 
in law or rather from a pathwise type convergence? 

The paper is organized as follows: the next section reviews some very recent results con- 
cerning the asymptotic behavior of weighted power variations of fractional Brownian motion. 
In section 3, after recalling the definition and the main properties of the so-called Newton- 
Cotes integral, we explain how to use it in order to study (jl.ip . Finally, in section 4, we state 
and prove our results concerning the exact rate of convergence associated to (11 . 131) . 
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2 Asymptotic behavior of weighted power variations 



Let k > 2 be an integer, and h, g : K. — * R be two functions belonging to C°°. Assume 
moreover that h and g are bounded with bounded derivatives. Denote by ^2n the 2re-moment 
of a random variable G ~ ^(0, 1). The following theorem collects some very recent results 
about the asymptotic behavior of the so-called weighted power variations of B, defined by 

71-1 

h(B e / n )(AB e / n ) K (recall that Afty n stands for B^+iy n - Bt/ n ). 

l=o 

Theorem 2.1 // k is even and H £ (0, 1) then, as n —> oo: 

n-1 „i 



£=0 

If K is odd and H E (0, !/2) then, as n — > oo; 

n— 1 



^DH-ij-^ ){AB r P»b _/Vfi / ^ (jBs)ds (2 2) 

5. // k is odd and H = 1/2 i/ien, as n -> 00, 

/ n-1 n-1 \ 

( B u n^ Y,9( B ^)( A V^ n ^ '£HB t/n )(AB t/n r+ 1 

e=o J 

B h J g(B s )(^dW s + fx K+l dB s ),fi K+1 J h(B s )ds), (2.3) 

with W another standard Brownian motion independent of B. 
[. If k is odd and H £ (V 2 i 1) then, as n — > 00: 

" _1 _ . /-Bi 



=0 £=0 
Law 



Y,HB e/n )(AB e/n r ^ M*+i / h{x)ds. (2.4) 

Remark 2.2 1. For sake of conciseness, we omit the proof of Theorem 12.11 We give below 
some ideas and references for the proofs. 



2. The convergence (12. ip is actually almost sure. Its proof is a classical result when h = 1 
(see e.g. [7] when k = 2). If /i is arbitrary, the proof could be completed along the lines 
of the proof of Lemma 3.1, p. 7-8 in [3]. 

3. Proofs of (|2.2p and (|2.4p can be completed along the lines of [17, Corollary 2]. 

4. We shall see that, for the standard Brownian motion case, in order to study our Mil- 
stein's type schemes one needs the behaviour of the triplet (|2.3|) and not only the 
behaviour of the second coordinate in (|2.3I) . The proof of (12. 3D can be completed along 
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the lines of [201 Corollary 2.9]. More precisely, using the methodology introduced in 
this latter reference, we first prove thalHI, as n — > oo, 



n-1 

K-l 



(Bt)te[o,i] > n ^ E g(Bi/n)(AB e/n ) K J 

(Bt)te[o,i], J q 9(Bs){VJ^dW s + » K+1 dB s )y (2.5) 



1=0 

Law 



Then, using the fact that (see (j2.1h ). as as n — > oo, 

n-1 „i 



and that Jq h(B s )ds is a random variable measurable with respect to B, the desired 
conclusion follows easily. 

5. Other results on weighted variations of fractional Brownian motion (or related processes) 
can be found in pQ and [9j. 

3 Newton-Cotes integral and fractional SDEs 

In the sequel, we will use, as integral with respect to B, the so-called Newton-Cotes integral 
introduced in [5] and studied further in [19] : 

Definition 3.1 Let f : R — > R be a continuous function, X,Y be two continuous processes 
on [0, 1] and N £ NU {0}. The N -order Newton-Cotes integral of f(Y) with respect to X is 
defined by: 

f f(Y s )d NC > N X s :=limprob- f ds{X s+s -X s ) f f(Y s + a(Y s+£ - Y s ))u N (da), t € [0,1], 
Jo £ i° e Jo Jo 

(3.1) 

provided the limit exists. Here vq = Sq, v\ = (*o+«5i)/2 and, for N > 2, 

2 ^T ( /• 1 TT 2(jV-l)n-^ \ 

1=0 \ J0 k*j J ~ k ) 

5 a being the Dirac measure at point a. 

Remark 3.2 1. The 0- and 1-order Newton-Cotes integrals are nothing but the forward 
integral and the symmetric integral in the sense of Russo-Vallois [23], respectively: 



f f(Y s )d NC >°X s = f f(Y s )d-X s = limprob- f f(Y s )(X s+£ - X s )ds, 
Jo Jo £ l° e Jo 



Using the notion of stable convergence for random variables, (|2.5j) is equivalent to say that 

n ^ £ 9(Bi/ n )(AB e/n ) K ^ B Z^ hl y f 1 g (B s ){^dWs + u K+1 dB s ). 
e=o Jo 

Here, denotes the cr-field generated by (B*)te[o,i] (see also Theorem 1.1, p. 3 in [6]). 
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and 



[ fiYJd^Xs = [ /(y s KX s =limprob- / f<Xa+e) + f(X '\ x a+e -X a )d8. 
Jo Jo £ i° £ Jo 2 

2. Another way to define vn is to view it as the unique discrete signed probability carried 
by j/(2A/-2) (j = 0, . . . , 2N — 2), which coincides with Lebesgue measure on polynomials 
of degree smaller than 2N — 1. 



The Newton-Cotes integral defined by (13. If) is actually a special case of so-called iV-order 
^-integrals introduced in [5], p. 789. Moreover, in the same cited paper, p. 795, one proves 
that the iV-order Newton-Cotes integral of f{B) with respect to B exists for any / G C 4N+1 
if and only if H G ( 1 /(4iV+2), 1). In this case, an Ito's type change of variables formula holds: 
for any antiderivative F of /, we can write 

F(B t )-F(0)= f f(B s )d NC ' N B s , tG [0,1]. (3.2) 
Jo 

Moreover, as a consequence of (|3.2p . let us note that 

/ f(B s )d NC ' N B s = I f(B s )d NC > n B s = F(B t ) - F(0), 
io Jo 

as soon as / G C AN+1 , n < N and H G ( 1 /(4n+2), 1). Therefore, for / regular enough, it is 
possible to define the Newton-Cotes integral without ambiguity by: 

f f(B s )d NC B s := / f(B s )d NC ' n B s if He (V(4„+2), 1). (3.3) 

Set n H := inf{n > 1 : H > 1 /(4n+2)}. Hence, an immediate consequence of (|3.2|) and (|3.3|) is 
that, for any H G (0, 1) and any / : R — > M. of class C 4n ^ +1 , the following Ito's type change 
of variables formula holds: 

F(B t ) = F(0) + [ f(B s )d NC B s , for any antiderivative F of /. (3.4) 
Jo 

Remark 3.3 In the sequel we will only use the fact that the Newton-Cotes integral verifies 
the classical change of variable formula (|1.3p or (13. 4p . Consequently, any other stochastic 
integral verifying (|3.4|) could be used in the following. 

All along this paper we will work with an ellipticity assumption, and we will also need 
regularity for the function a. More precisely, we suppose 

((f) inf]g | cr | > and a G C°°(R) is bounded with bounded derivatives. 

Under hypothesis (<f), the flow associated to a, given by (|1.2|) . is well-defined and verifies 
the group property: 

\/x,y,zeR, 4>(cf)(x,y),z) = 4>(x,y + z). (3.5) 
Note that the process X x given by (jl.4p verifies: 

Xf = x+ I a(X^)d NC B s , t G [0, 1], (3.6) 
Jo 

as we can see immediately, by applying (13.41) . 
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Remark 3.4 In |19| (see also |21|). one studies a notion of solution for (|3.6p and also the 
existence and the uniqueness of solution. Note however that, in the present work, we will 
only use the fact that there exists a natural solution to (13, 6p given by (jl.4j) . 

The following result explains the definition (j!.13[) . By using (j 1 . 1 2 f) , the process X x defined 
by (|1.4p can be expanded as follows: 

Lemma 3.5 For any integers m > 0, n > 1 and i G {0, . . . , n — 1}, we have 



3=0 U 



3+1 



+ d™B tl d™B t2 ... d™B tm+1 ® m+l «{X? m+2 )d^B tm+2 . 

Ji/n Jt/n Jt/n Jt/n 



(3.7) 



Proof. We proceed by induction on m. By applying (|3.6|) . and then using (jl,4|) and (|3.4|) . 
we can write: 



X? w)/n = Xf /n + a(Xf /n )AB e/n + / {a{XD - a(Xf /n ))d* c B h 

= Xf /n + a(Xf /n )AB e/n + / d" c B tl / aa\Xl )d NC B t , 2 , 



which is exactly (|3.T[) for m = 0. 

Now, let us assume that (13, 7p is true for some m 6 N U {0}. Then we can write 

m 1 



+ ^ m+ V(X^J / d NC S 4l . . . / d" c B tm+1 I d" c B tm , 2 



+ / d^B tl ...j d NC B tm+1 j {@™+MXL +2 )-® m+1 °(X? /n ))d" C B tm+2 . 



(*+l)/n ftm ftm+l 

f d NC B tm+1 

t/n Jl/n 

(3.8) 

On one hand, using (|3.4p repeatedly, it is immediate to compute that 

r(t+l)/n ft m ft m +l i , _ 

/ d^B tl ... d^B tm+1 d™B tm+2 = (AB e/n r+\ 

Jl/n Ji/n Ji/n [m + Z)\ 

On the other hand, using (jl.4p and again ()3.4p . we can write 

rtm+2 



rtm+2 

® m+1 a(X? m+2 ) - @ m+1 a(Xf/J = / a{9 m+l a)\xt m+3 )d^ c B 

J £/n 



m + 3 ' 



Finally, putting these latter two equalities in (13. 8p and noting that o-(2$ m+1 a)' = & m+2 a by 
definition, we obtain that (|3.7|) is true for m + 1. The proof by induction is done. □ 

Clearly, (|1.13p is the natural scheme constructed from (|3,7p . by considering the (m + 2)th 
multiple integral in the right hand side of (|3.7f) as a remainder. 



7 



4 Rate of convergence of the approximating schemes 
4.1 Statement of the main result 

Recall that we denote by fi2n the 2n-moment of a random variable G ~ ^(0, 1). For m £ 
N U {0}, let us introduce the functions g m , h m : R —* R given by: 



5// 



-a' h m + h m+ i and /i r , 



(^ ra+1 cr)/ (7 

(m + 2)! ' 



(4.1) 



Our main result contains a complete answer to the questions in the introduction and can 
be stated as follows: 

Theorem 4.1 Assume that hypothesis {£) is in order, and let m £ N U {0}. Then, for any 
H £ ( 1 /{m+2),l), the sequence {X^ }„ 6 n defined by tl.l3\) converges almost surely toward 
Xf = 4>(x, B\) as n — > oo. Moreover, 

• when m is even and H £ (l/(m+2), 1), 



n 



(m+2)H-l 



Prob 



fi m+2 a(Xf) C h m {Xf)ds- (4.2) 
J o 



when m is odd and H £ ( 1 /(m+2), !/2), 



(m+3)ff-l 



x{ n) - Xf 



Prob 



-cr 



^)(Xf)d6 



u>/ien m is odd and H = l /2, 



n 



(m+l)/ 2 



X 



(«) 



Xf 



Law 



ff(Xf ) / h m {X X s )y ^2m+A dW s + ^ m+3 d5 s 



+|"m+3 / 9m(Xs)ds 



with W a Brownian motion independent of B. 
when m is odd and H £ (!/2, 1), 



■n 



(m+l)H 



x[ n) - Xf 



Prob 



fJ-m+3 O-(Xf) 



Bi 



hm{4>{x,y))dy. 



Remark 4.2 1. For m = and -ff > 1/2, one recovers the convergence (11.6 



(4.3) 



(4.4) 



(4.5) 



2. With the same method used to obtain (14. 4p . one could prove (II. lip with the help of 
Lemma 3.4 in |19j . Details are left to the reader. 

3. Actually, we could prove that the convergence is almost sure in (|4.2[) . Also the conver- 
gence in (14. 3p is certainly almost sure, but the method of proof we have used here does 
not allow to deduce it. Thus it remains an open question. 

4. According to Theorem l4.ll whenever H £ ( 1 /(m+2), 1 /2) the scheme X of size m = 2k — 1 
has the same rate of convergence than the scheme X of size m = 2k, namely nP ,K+ ' i ' H ~^ . 
Thus, it is a priori better to use odd-size schemes. 
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5. With the same method (see also Theorems 2 and 4 in [15]). one could also derive the 
exact rate of convergence for the global error on the whole interval [0, 1]. 



Observe that, under the convergences (|4.2h and (|4.3p give the right lower error bound 
if the probability that the right-hand side vanishes is strictly less than 1. Due to (|1.4|) and 
the fact that Bt has a Gaussian density for any t €]0, 1], it is easy to see that this last fact 
is equivalent to have that the real function inside the integral, say f m , is not identicaly zero. 
Indeed, if g(B s )ds = almost surely for a certain g G C£(R), then = D u g{B s )ds = 
g'(B s )ds, for any u G [0,1] (here D denotes the Malliavin derivative with respect to B). 
We deduce that g'(B u ) = 0, for any u E [0, 1], and, since the support of the law of B\ (for 
instance) is R, we obtain g' = 0. The desired conclusion follows easily. 

Except for m = 0, solving f m = seems complicated. Nevertheless, when m = 1, we can 
state: 

Proposition 4.3 Assume that {$) is in order, and moreover that a does not vanish. Then 
the function 3cr' 3 + 6aa'a" + a 2 a'" (which is, up to a constant, the function appearing inside 
the integral of the right-hand side of |^.3| ) when m = 1) is not identically zero. 



Remark 4.4 1. When o~(x) = a is constant, we have x[ n ^ = Xf = x + a B\. Conse- 
quently, the study of the rate of convergence in the case where a is a constant function 
is not interesting. 

2. A corollary of Theorem 14 . 1 1 and Proposition 14.31 is that, under the additional hypothesis 
that a' does not vanish, the upper and lower error bounds always come from a conver- 
gence in probability whenever H ^ l /2 and m = 1. In particular, we never observe a 
phenomenon of the type (jl.9p . 



Proof of Proposition 14.31 Since a' does not vanish, we have either a' > or a' < 0. 

Suppose for instance that a' > 0, the proof for the other situation being similar. Assume for a 
moment that / := 3c/ 3 + 6a a' 'a" + a 2 a'" is identicaly zero. We then have 3a' (a 2 )" = -a{a 2 )"' . 
We deduce that the derivative of a 3 (a 2 )" is zero and then (a 2 )" = a a~ 3 on R, for some 
Set h = a 2 ; we have h" h' = ah' /i -3 / 2 or, equivalently, h' 2 = —4a /i -1 / 2 + (3 for some /3 & R. 
In particular, we have (3 — 4a y~^ 2 > 0, for any y £ /i(R). Let F be defined on /i(R) by 



F{y) 



dz 



lh(o) v 7 /? - 4a z- 1 / 2 
For all x £ R, we have 

F(a(x) 2 ) = F(h(x)) = x + 7, for some 7 G R. (4.6) 

The function a 2 being bounded, we necessarily have h(x) — > ( 4a //3) 2 (in particular [3 7^ 0) as 
x — > 00. Then, since h" = ah~ 3 ^ 2 , this implies that h"(x) — > P 3 /(a 3 oi 2 ) as x — > 00, which is in 
contradiction with the fact that h = a 2 is bounded. The proof of the proposition is done. □ 



4.2 Proof of Theorem 14.11 

Here, and for the rest of the paper, we assume that H belongs to ( 1 /(m+2), 1) and we denote 
A n = maxfc^.^n-i \AB k / n \. We split the proof of Theorem 14.11 into several steps. 
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1. General computations. The following lemma can be shown by using the same method 
as in the proof of Lemma 13.51 but with Lebesgue integral instead of Newton-Cotes integral 
(and by taking into account that a G C°°(M) is bounded with bounded derivatives, in order 
to have uniform estimates): 

Lemma 4.5 As y — > ; we have, uniformly in x 6K, 

m+2 

( / ) (x,y)=x+Yl 7— — jT7 @ j o-(x) y j+1 + 0{y m+A ). 



By applying this lemma to x = X k ^ and y = AB k / n , we obtain, using the definition of 

(fc + l)/n' 



A (fe+1 



X. 



(n) 

(*+l)/n 



^xijl,AB k/n )-^a(X i 



(AA /n ) 



m+2 



(m + 2)! 



t+ 2..v(^( A %) W+3 



/" (m + 3)! 



+o(Ar 4 )- 



By straightforward computations we getH 

= ^ > + h m (X$)(AB k/n ) m+2 + fa (^)(A ft/ „) m+3 + Q(K 



(4.7) 



with g m and h m given by (|4.1|) . By applying the group property (|3,5|) repeatedly, we finally 
obtain that, for any I E {1, . . . , re}: 



X$ = <j> x ,B Vn + ^ h m (X$)(AB k/n ) m+2 + J29m(xljl)(AB k/n ) m+3 + 0(nA 

V fc=0 fc=0 

Since 9<t>/dy = a o is bounded and (<f) is in order, we deduce, as re — > 00, 



sup 

ee{i,...,n} 



y(n) yx 



sup 

ee{l,...,n) 



X 



(n) 



0(nA™+ 2 ). 



(4.8) 
(4.9) 



In particular, converges almost surely to Xf as re — » 00, since -ff > 1 /(m+2). 



Proof of fiJJty . Let m be an even integer. As a consequence of (I4.8|) and (14. 9p . we can 

write 



X 



(n) 



U B X + J] M*yJ (Ai? fe / n ) m+2 + 0(nA-+ 3 ) + 0(n 2 A^+ 4 ) J . (4.10) 



In fact, we rather obtain 



X%i )/n = (x$,AB k/n + h m (xijl)(AB k/n ) m+2 



•S»(^)(ABv 



•o(a:: 



which is not exactly (|4.7[) . But, in order to introduce 0(A™ +4 (_B)) in the argument of <j>, we proceed as follows, 
by using the ellipticity property in hypothesis (S ): 

tj>(x, z) + 0{5) = <j>(x, (j>~\x, 4>{x, z) + 0(5))) = <j>(x, z + 0(S)). 
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(Observe that the main difference between (|4.8p and the previous identity is that the argument 
of h m is here X®,). By using (|2.ip with k = m + 2, and due to the fact that Xf = (f>(x, B t ) 
and d <f>/dy = & o (j)^ we finally obtain (j4.2jl . 



5. Proof of \4-3\ ) for H G ( 2 /(2m+3), !/2). Let m be an odd integer and assume that 
H G ( 2 /(2m+3), V 2 )- Thanks to (|4.9|) . identity (|4.8|) can be transformed into 



4"j = ( x ,B Vn + Y j h m (X^)(AB k/n ) m+2 

k=0 

+ ^ 5m (^ / J(A^) m+a + 0(nAr 4 ) + 0(n 2 A^ 5 ) ) . (4.11) 



j m+3 

k=0 / 

On the other hand, due to (<f), we have, for any fixed M > 1 and uniformly in x E 

Ai 



^(x, y2 ) = </>(*, yi ) +Y,-M( x > - + o((»2 - yi) M+1 ). 

Combined with (|4.8p . it yields 

= + £ ? ( E m*#j ( A ^/n) m+2 + o(nAr 3 : 

3=1 J ' y \ki=0 

+ 0(n M+1 A( t m+2 )( M+1 )). (4.12) 
By using (|4.12p with M = 1 as well as the equality d^/dy = a o 0, we get 

= + *(*?/J E ^(X^JCAS^)^ 2 + 0(n 2 A^ 4 ) + 0(nAr 3 ) 

fci=0 

and then, by (fOI) . 

= *V» + E ^(^ l/?i )(Ai? fcl/ri ) m+2 + 0(n 2 A 2 ™+ 4 ) + 0(nAr 3 )- 

fci=0 

By inserting the previous equality in (|4.11|) with t = n, we obtain 

n—l n—1 

X[ n) = <j>{x,B 1 + £ h m (XZ /n ){AB k/n ) m+2 + Y,9m(X- /n ){AB k/n ) m+3 

k=0 k=0 
n-1 fc-1 

+ E (^r +2 E M*7/J (^/X +2 

fc=0 j=0 

+ 0(n 3 A3 m + 6 ) + 0(n 2 A 2m+5 ) + 0(nA™+ 4 )). (4.13) 
Due to (12. 2|) with k = m + 2 we have, as n — > 00, 
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and also, due this time to (|2.ip with k = m + 3, as n ^ oo, 

n-l „i 

n (™ + 3)H-iJ2 9m (X* /n )(AB h/n ) m+3 ^ b Mm+3 / g m (X*)ds. 
k=0 Jo 



Moreover, since we assume in this step that H > 2 /(2m+3), we have, as n — > oo, 

i(m+ 3)H A rn + 4 P£ob Q) n ( m+3 )// + l A 2m+5 Prob Q ^ n ( m+ 3)H+2 A 3 m+ 6 Q 



At this level, we need the following result which is contained in |17[ Proposition 7]: 

Lemma 4.6 Fix an integer q > 2 and denote by H q the qth Hermite polynomial. Let f £ 
C 2q (M.) be bounded with bounded derivatives and, for k £ {!,••• ,n}, denote 



k-l 



3=0 



Then 



E\S^\f)( 2 = 0(n lv ( 2 2Hq ^ as n —>■ oo, uniformly in k. 



(4.14) 



Recall also that, since m + 2 is odd, the monomial x m+2 may be expanded in terms of the 
Hermite polynomials as follows: 



m+2 _ 

H"in-\ 



(m + l)/ 2 

l+3 x = E a m+ 2,2q+iH2q+i(x), for some universal constants a m+2 ^ q +i ■ 

9=1 

(4.15) 

Therefore, for k G {1, . . . , n}, 



k-l 



11 



3=0 



m - fJ-m+3 n 2H ~ 1 '^2h m {Xff n ) ABj/ n 

3=0 



(m+l)/ 2 

n H ~ l E a m+2 ,2g+iS^ q+1 \h m (cl)(x,-))) 

q=l 



and, by KW\ . 



E 



k-l 



k-l 



n 



(m^ 



+3)H - 1 E M*?/J (A^/,) m+ " n 2H - 1 E h m (Xf /n )AB. 



i/n 



Hence 



E 



3=0 



n—1 



3=0 



( n <2H-lM-4H)y 



k-l 



,m+2 



fc=0 



i=o 



n— 1 



fc-1 



A:=0 



3=0 



^(|-fl-mI)V(l-4I-mff)j 
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which tends to zero as n — ► oo, because H > 2 /(2m+3) implying £f > 1 /(2m+2) and H > i/(m+4). 
Moreover, by the mean value theorem: 



fc-i 



f Bk/n h m {4>{x,z))dz - ~J2<rh' m (Xl /n )(AB. 

J ° 3=0 



j/n) > 



for some 9j/ n between i/n and (i+ l )/n. Consequently, since H < 1/2, we have 



E 



k-l 

J2h m (Xf /n )AB j/r 

3=0 



0(n 



2-4#n 



so that 



E 



n-1 k-l 
,2if-l^.i.( Ivx \(\r> \m+2 



n 



1 Y, °ti m {XZ /n ) (AB k/n ) m+ Y, h m {X* /n )AB j/r , 

3=0 



k=0 



0(n 



-(m+2)H\ 



0. 



Finally, by combining all these convergences to zero together, we get 



7V 



n-1 k-l 
v m+2 ' 



m+2 Prob 



k=0 j=0 

so that the proof of (|4.3p is done in the case where if > 2 /(2m+3). 

4- Proof of (f^.ffl ) /or if G (Y(m+2), 2 /(2m+3)]. It suffices to use (14. 12f) with the appropriate 
M for the considered H and then to proceed as in the previous step. The remaining details 
are left to the reader. 

5. Proof of By going one step further in (|4.1ip using (|4.9p . we get 



n-1 



n-1 



X[ n) = </>(x,B 1 + Y J h m {X x k/n ){AB k/n ) m+2 + Y J 9m{X x k/n ){AB k/n 



, m+3 



k=0 



k=0 



+ 0(n 2 A 2 n m+4 ) + 0(nA™+ 4 )). 



Whenever m > 3 and since H = 1 /2, we have, as n — > oo: 



(m + 



D/2+1^+4 Prob Q ^ n ( m+1 )/ 2+ 2 A 2 m+ 4 ^ob Q _ 



Hence, for m > 3, (|4.4p is an immediate consequence of (|2.3p and of the previous two relations. 
If m = 1, we rather need to use (|4.13p . Since H = 1/2, we have, as n — > 00: 

n 2 A 5 P£ob Qj n 3 A 7 Q ^ ^9 ^ob Q 

Finally, combining these convergences with (12. ip (for ii = 1/2), (12. 3p and the fact that 



E 



n-1 



fc-1 



n 



Y ^1 {X x k/n ) (AB k/n y £ hi (Xf /n ) (AB 



j/n 



k=0 



3=0 



n-1 k-l 



n a 2X)EaJii(A' Vs )(%) d /i 1 (A? / J(AB j/ ,) i = 0( 

fc=0 j=0 



7? 



as n — > 00, 
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we obtain (|4.4p also for m = 1. 

6. Proo/ 0/ {J3p. By combining (|4.10[) with the fact that Xf = 4>(x, B t ) and 9<t>/dy = croft, 
we get 

^ = *? + a(Xf)J2h m (X^ /n )(AB k/n ) m+2 + 0(nA^) +0(n 2 A^+ 4 ). 

fc=0 

Since H > 1 /2, we have, as n — > 00, 

n(m+1) H + l A ™+3 P^b Q ^ n(m+ l)H+2 A 2 m+ 4 Q 

Hence (I4.5P is an immediate consequence of (I2.4p . _ 

Acknowledgement: We thank Jean Jacod for helpful discussion about stable convergence 
and especially about reference [6] and Serge Cohen for its suggestion of reference [7J. 
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